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Precision measurements of the polarization of the cosmic microwave background (CMB) radiation, 
especially experiments seeking to detect the odd-parity "B-modes", have far-reaching implications 
for cosmology. To detect the B-modes generated during inflation the flux response and polarization 
angle of these experiments must be calibrated to exquisite precision. While suitable flux calibration 
sources abound, polarization angle calibrators are deficient in many respects. Man-made polarized 
sources are often not located in the antenna's far-field, have spectral properties that are radically 
different from the CMB's, are cumbersome to implement and may be inherently unstable over 
the (long) duration these searches require to detect the faint signature of the inflationary epoch. 
Astrophysical sources suffer from time, frequency and spatial variability, are not visible from all 
CMB observatories, and none are understood with sufficient accuracy to calibrate future CMB 
polarimeters seeking to probe inflationary energy scales of 10 GeV. Furthermore, both man-made 
and astrophysical sources are often much brighter than the CMB B-mode signal and these bright 
sources can cause non-linearities in the detector's response. Both man-made and astrophysical 
sources require dedicated observations which detract from the amount of integration time usable for 
detection of the inflationary B-modes. CMB TB and EB modes, expected to identically vanish in the 
standard cosmological model, can be used to calibrate CMB polarimeters. By enforcing the observed 
EB and TB power spectra to be consistent with zero, CMB polarimeters can be calibrated to levels 
not possible with man-made or astrophysical sources. All of this can be accomplished without any 
loss of observing time using a calibration source which is spectrally identical to the CMB B-modes. 
The calibration procedure outlined here can be used for any CMB polarimeter. 
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Introduction: Inflation is perhaps the most promis- 
ing model of the early universe, resolving the Big 
Bang model's flatness and horizon problems and pro- 
viding seed perturbations for structure formation (see 
e.g., [JJ for review). Besides the density, or "scalar", 
seed perturbations, inflationary cosmological models 
also predict "tensor" perturbations arising from a 
primordial gravitational wave background. Primor- 
dial scalar perturbations create only CMB B-modes, 
while primordial tensor perturbations generate both 
parity-even B-modes and parity-odd .B-modes polar- 
ization [2H1]. The detection of primordial tensor B- 
modes in the CMB would confirm the existence of 
gravitational wave perturbations in the early universe. 
Numerous observational efforts are underway to de- 
tect the CMB's -B-mode since such a detection would 
establish the energy scale at which inflation occurred. 
The amplitude of primordial .B-modes can be char- 
acterized by the tensor-to-scalar ratio, r. The most 
restrictive limit on r is currently r < 0.18 (95% confi- 
dence) [5] , and the best direct limits on r from B-mode 
measurements is r < 0.72 (95% confidence) [5J. 

An impressive battery of CMB experiments have 
mapped the CMB's intensity to near cosmic- variance- 
limited precision to I ~ 3000[5J [7J |S]. However, the 



measurement of the CMB's B-mode polarization sig- 
nal is considerably more challenging since it is 10 
to 100 times smaller than the CMB's temperature 
anisotropy. Compared to these signals the .B-mode 
polarization from inflationary gravitational waves is 
even more challenging - current upper limits corre- 
spond to B-mode fluctuations less than 10% of the 
E-modes. Calibration accuracy required by a given 
instrument to constrain or detect the minute B-mode 
signal is well-discussed in the literature [SJ. Calibra- 
tion is accomplished either using hardware calibrators, 
located in the near field of the instrument [PTjHT5] , or 
by measurements of polarized astrophysical sources, 
e.g. the Moon, Tau A, Cen A, 3C 273, or the galactic 
plane [Mlfi"7] . Neither hardware polarization calibra- 
tors nor astrophysical sources can achieve better than 
~ 1° precision on the polarization angle calibration. 
It is hard, if not impossible, to do better, yet a pre- 
cision of 1° is insufficient for detecting r = 0.025 to 
0.01 [9] - the goal of future polarimeters. Rather, it 
has been shown [TB] that for r = 0.01 to be biased 
by < O.loy (where oy is the nominal statistical uncer- 
tainty on the inferred value of r), the uncertainty in 
pixel rotation cannot exceed 4' even for a nearly ideal 
CMB experiment. Even relaxing this requirement, by 



allowing a bias on r as large as ~ lay, only permits a 
miscalibration level of < 12'. 

Miscalibration of the instrument's polarization an- 
gle (pixel rotation) mixes polarization modes, leaking 
E- into B-modes, thereby producing spurious B-mode 
polarization [§J H3 [2U] . Additionally, due to polariza- 
tion mode-mixing, new TB and EB correlations are 
generated. Since the standard cosmological model is 
parity-even the TB and EB correlations identically 
vanish. Therefore the TB and EB spectra can be used 
to probe the miscalibration of the pixel rotation an- 
gle. Furthermore, the miscalibration angle itself can 
be quantified, then removed, resulting in an unbiased 
measurement of B-mode polarization. This calibra- 
tion procedure is accomplished during data acquisi- 
tion, requiring no additional observing time. There- 
fore it is referred to as "self-calibration". Moreover, 
since the calibration signal is the CMB itself, any con- 
cerns that the detector response will behave nonlin- 
early are eliminated. 

It has been shown that miscalibration produces a 
distinct signature in the (otherwise zero) (EB) and 
(TB) correlations. The amplitude of these "forbid- 
den" correlations is proportional to the amount of mis- 
calibration and, furthermore, it is known that several 
other instrumental systematics can be detected using 
these EB and TB correlations [Hl[2Qj. This paper 
uses these correlations to calibrate CMB polarimeters 
to levels not achievable with laboratory or astrophys- 
ical sources. 

Polarization Map Making and Miscalibration: Fol- 
lowing [6], the timestream data from a single polar- 
ization sensitive detector, d,, is written as 

di = 9i [T(p) + 7i(Q(p) cos 2^ + U(p) sin 2^)] , (1) 

where g t is the flux calibration, or "gain" for the i- 
th detector, T, Q, U are the beam-integrated CMB 
Stokes parameters for the map pixel in direction p, 
7i = (1 — £j)/(l + e,) is the polarization efficiency fac- 
tor, e, is the polarization leakage for the i-th detector, 
and tpi is the detector's polarization orientation pro- 
jected on the sky. The goal of mapmaking is to recover 
T,Q,U from the detector timestreams. 
The angle ipi is modeled as 

Ipi = V'dosign + Alp , (2) 

where V'dosign is the intended orientation of the de- 
tector on the sky with respect to right ascension and 
declination and Aip is the miscalibration of the detec- 
tor. Figure 1 displays the coordinate system used for 
a single polarization sensitive detector. 

Calibration of the detector's pixel rotation is one 
of the most challenging tasks facing the experimen- 
talists [TTJ [3T] . A miscalibration of the detector's po- 
larization angle by an amount Aip rotates primordial 




FIG. 1: Coordinate system showing the relevant angles 
from equation (2) for a single polarization sensitive detec- 
tor. The horizontal and vertical axes are right accession 
(a) and declination (5). For this detector i/'design was in- 
tended to be parallel to a. 

Stokes parameters Q(n), and U(n) into the observed 
quantities: 

Q(n) ± iU(n) = e ±2lA ^(Q(n) ± iU(n)) . (3) 

The measured Fourier modes, E(l) and B(l), written 
in terms of the primordial modes E(l) and B(l), be- 
come 

E(l) = cos(2A^)£(l) + sin(2AV0-B(l) 

B(l) = - sin (2Aip)E(l) + cos (2At/>)B(1) . (4) 

The above equations show that pixel rotation mod- 
ifies the power spectra of E(l) and B(l) and generates 
spurious correlation between E(l) and B(l) and be- 
tween T(l) and B(l), modifying the observed power 
spectra as follows: 

Cj E = cos (2AiP)Cj E 

Cf E = sin 2 (2Aip)C BB + cos 2 (2Aip)Cf E 

Cf» = l r \n(^)(C? B -C EE ) 

Cj B = - sin (2AiP)Cj E 

Cf B = cos 2 {2AiP)Cf B + sin 2 (2AiP)Cf E . (5) 

Here, and throughout, tildes represent primordial 
quantities. From Eq. [5] it is clear that pixel rota- 
tion generates spurious -B-modes in the absence of 
primordial B-modes. Calibration of the pixel rotation 
involves finding Aip for the detector system and re- 
moving it from the data prior to map making (Eq. [T]) . 
This procedure will recover the unrotated, primordial 
CMB polarization spectra. 
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FIG. 2: Simulation demonstrating the self-calibration of pixel rotation using the CMB's EB power spectrum. The 
polarization angle's miscalibration angle is assumed to be Atp = 2°, and the instrumental noise is Ae — Ab = 5[iK- 
arcmin. The beam's full- width at half-maximum is Bfwhm = 4'. The solid red curve shows the _E-mode power spectrum, 
the green long-dashed curve shows the B-modes induced by miscalibration. The blue short-dashed curve shows the B- 
mode after correcting for the miscalibrated pixel rotation angle. Finally, to demonstrate how well the self-calibration 
works, the instrumental noise contribution is subtracted from the de-rotated £?-mode in the bottom-most (magenta 
dotted) curve; these are residual B-modes remaining after self-calibrating the instrument. The primordial B-mode 
spectrum corresponding to r = 0.01 is shown in dot-dashed (cyan) line. For clarity, the lensing B-modes are omitted. 



In the absence of cosmological parity violation lead- 
ing to cosmic birefringence [351 123] > the detection of 
EB and TB spectra - each having the same implied 
pixel rotation angle Aif) ~ directly implies that the de- 
tector's polarization angles have been miscalibrated. 
Such miscalibration can be caused, e.g., by fabrication 
errors or sources of birefringence in the telescope's op- 
tics. 

The miscalibration angle Atp, is obtained from the 
observed Cj B and Cf B by minimizing the variance 
of the difference between the observed and theoretical 
power power spectra as a function of A.ip. 

Using Eq. [5j the two independent likelihood func- 
tions for the miscalibrated pixel rotation angle, Aip, 
become 



For simplicity, in Eq. [fj] it is assumed that Cf E 
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where ] sky is the fraction of the sky observed, 
{T, E, B}, and 
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Here, Ax is the detector noise and Ofwhm is the full- 
width at half-maximum (FWHM) resolution of the 
polarimeter's Gaussian beam. 

The best-fit estimates for the pixel rotation angle 
are obtained by maximizing the likelihood functions, 
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Eq. [6j resulting in 

Alp TB = 
AlpEB = 
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and x e {t,^,^}. 

The values for AlpEB and AlpTB obtained from ob- 
servations, Eqns. [9j should agree to within the statis- 
tical uncertainty (which is derived next). The consis- 
tency AlpTB — AlpEB provides a powerful cross-check 
on the hypothesis that the pixel rotation has been 
miscalibrated. Since there is intrinsic correlation be- 
tween EB and TB, for real data Eqns. [6] should be 
replaced with a single likelihood function that simul- 
taneously uses C TB , C EB and their correlation, to 
infer a single Aip value. In practice, AlpEB will be su- 
perior to AtpT b but the value of having two indepen- 
dent estimates, which must agree in order to apply the 
self-calibration method, motivates the construction of 
both estimators. 

The signal-to-noise for the EB and TB detection is 
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where f s k y is the fraction of the sky observed. The 
uncertainties in AlpEB and AlpTB are 
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Equations [9] and [12] give the miscalibration of pixel 
rotation and its uncertainty. Armed with these quan- 
tities, the experimentalist can go back and correct the 
assumed polarization angle by subtracting Aip from 



ipi in Eq. ^ For a toy experiment with A p — 5/jiK- 
arcmin, f s k y = 0.01, and Qfwhm = 4' a pixel rota- 
tion as small as Aip ~ 0.05° (3') can be detected using 
the EB power spectrum alone, and Aip ~ 0.1°, using 
the TB power spectrum alone. These values are more 
than sufficient to detect r = 0.01 with < O.loy bias. 
Fig. [2] shows the pixel rotation calibration method de- 
rived using the EB spectrum. A similar result can be 
obtained for Atp TB - 

Discussion: While pixel rotation is potentially the 
most pernicious obstacle to detecting primordial B- 
modes, there are other systematic effects such as dif- 
ferential ellipticity and differential pointing that can 
produce TB and EB correlations [HI GO] • However 
the TB and EB correlations induced by these system- 
atic effects have different angular dependencies from 
pixel rotation [24 . Therefore, high-sensitivity obser- 
vations covering large numbers of multipoles can make 
detailed measurements of the TB and EB correla- 
tions to quantify, and correct for these types of sys- 
tematic errors. Furthermore, due to the dependence 
on the beamsize, for a fixed level of instrument noise, 
a higher resolution experiment will always calibrate 
pixel rotation more precisely than a lower resolution 
experiment. 

Pixel rotation and cosmological birefringence are 
fully degenerate effects; i.e. the angle Aip in Eq. ^ 
could be the sum of the two effects. However, for the 
purpose of B-mode detection it does not matter what 
causes the polarization rotation Aip; we simply de- 
rotate the polarization by the angle Aip inferred from 
the EB and TB estimators. Furthermore cosmic bire- 
fringence detection is not sacrificed when cosmic bire- 
fringence is spatially varying [20l I24H28] . 

Complications arising from the E — B separation 
due to partial sky coverage and the EB and TB 
correlations induced by gravitational lensing of the 
CMB by large scale structure are irrelevant to the 
self-calibration procedure proposed here because these 
EB and TB correlations couple only different ^-values, 
i.e., (EeBgi) and [TiB^) will be non-vanishing only 
when i^£'. 
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